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ABSTRACT 

The  foundation  of  the  theory  of  functionally  graded  plates  with  all  four  edges  simply  supported, 
under  a  Eriedlander  explosive  spherical  air-blast,  is  developed,  within  the  classical  plate  theory 
(CPT).  The  constituent  materials,  ceramic  and  metal,  vary  across  the  wall  thickness  according  to 
a  prescribed  power  law.  The  theory  incorporates  the  geometrical  nonlinearities,  the  dynamic 
effects,  compressive/tensile  edge  loadings,  damping  effects,  and  the  structural  symmetries 
(symmetric  and  asymmetric).  The  static  and  dynamic  solutions  are  developed  leveraging  the  use 
of  a  stress  potential  coupled  with  the  Extended-Galerkin  method  and  the  Runge-Kutta  method. 
Validations  with  simpler  cases  existing  within  the  literature  are  made.  The  analysis  focuses  on 
how  to  alleviate  the  unwanted  effects  of  large  elastic  stresses  and  deformations  through  material 
selection  and  proper  gradation  of  the  constitutive  phases. 

Key  Words:  Eunctionally  Graded;  Explosive  Blast,  Dynamic  Response,  Elastic  Stress 

1.  INTRODUCTION 

During  military  conflicts,  the  structure  of  army  military  vehicles  may  have  to  withstand  blast 
loading,  the  residual  effects  of  which  are  the  blast  pressure,  fragmentation,  and  heat.  Advances  in 
functionally  graded  materials,  which  combine  the  properties  of  two  dissimilar  materials,  has  been 
a  motivating  factor  in  viewing  these  types  of  materials  has  a  viable  alternative  to  the  current 
metallic  structures  being  utilized.  Eunctionally  graded  materials  (EGM’s)  are  high-performance 
heat  resistant  materials  able  to  withstand  ultrahigh  temperatures  and  extremely  large  thermal 
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gradients.  The  FGM’s  are  mieroscopically  inhomogeneous  with  thermo-meehanieal  properties 
whieh  vary  smoothly  and  eontinuously  from  one  surfaee  to  another.  These  graded  structures 
allow  the  integration  of  dissimilar  materials,  like  ceramics  and  metals  that  combine  different  or 
even  incompatible  properties  such  as  hardness  and  toughness.  In  addition,  by  varying  both  the 
material  selection  and  the  prescribed  variation  of  the  two  constituent  materials  from  one  surface 
to  another,  it  is  believed  that  the  excessive  stresses  and  deflections  can  be  diminished  resulting  in 
many  benefits  to  the  survivability  of  the  soldier. 

These  facts  necessitates  the  need  to  gain  an  understanding  of  the  response  behavior  and  the 
load  carrying  capacity  of  plate-type  structures  composed  of  both  symmetric  and  asymmetric-type 
symmetries  of  functionally  graded  plate -type  structures  exposed  to  blast  loading  and  to  gain  an 
understanding  of  how  varying  both  the  material  selection  and  the  prescribed  variation  of  the  two 
constituent  materials  from  one  surface  to  another  can  alleviate  the  excessive  mechanical  stresses 
and  deflections  occurring  during  a  blast  event. 

2.  BASIC  ASSUMPTIONS  AND  PRELIMINARIES 

The  plate  mid  surface  is  referred  to  a  cartesian  orthogonal  system  of  coordinates  {x,y,z),  while  z 
is  the  thickness  coordinate  measured  positive  in  the  upwards  direction  from  the  mid-surface  of 
the  plate  with  h  being  the  uniform  plate  thickness  of  the  plate,  and  y  is  directed  perpendicular  to 
the  x-axis  in  the  plane  of  the  plate.  See  Fig.  1  below. 

The  nonlinear  elastic  theory  of  FG  Plates  is  developed  using  the  classical  plate  deformation 
Theory  [6].  It  is  also  assumed  that  the  FG  plate  is  made-up  of  ceramic  and  metal  phases  whose 
material  properties  vary  smoothly  and  continuously  across  the  wall  thickness.  By  applying  the 
rule  of  mixtures,  the  material  properties  such  as  Young’s  Modulus,  Density,  and  Poisson’s  Ratio 
are  assumed  to  vary  across  the  wall  thickness  as 
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P{z)^P,V,{z)  +  PJ^^{z),  (1) 

in  which  Pc  and  Pm  denote  the  material  properties  of  the  eeramie  and  metallie  phases,  of  the 
plate,  respeetively.  F^(z)and  V^{z)  are  eorrespondingly,  the  volume  fraetions  of  the  eeramie 
and  metal  respeetively,  fulfilling  the  relation 


F,(z)  +  F^(z)  =  l.  (2) 

By  virtue  of  (2),  Eq.  (1)  ean  be  expressed  as 

P(z)^(Pc-PmWc(^)  +  Pm-  (3) 

By  observation,  one  ean  deduee  that  for  (z)  =  0 ,  P(z)  -  and  for  =  1 ,  P(z)  -  P^,.  As  a 


result,  E^(z)  e  [0, 1]. 

Two  Seenarios  of  the  grading  of  the  two  basie  eomponent  phases,  eeramie  and  metal,  through 
the  wall  thickness  are  eonsidered. 

Case  (1):  The  phases  vary  symmetrieally  through  the  wall  thiekness,  in  the  sense  of  having  full 
eeramie  at  the  outer  surfaees  of  the  plate  and  tending  toward  full  metal  at  the  mid-surfaee.  For 
this  ease,  V^.  (z)  ean  be  expressed  as 


Vc(z)  = 


yhllj 


l  +  sgp(z)^ 

V  2  J 


yhllj 


l-sgp(z) 

V  2  y 


where  the  Signum  funetion  is  defined  as 


sgn(z) = 


1,  z  >  0 
0,  z  =  0  . 
-1,  z  <  0 


(4) 


(5) 
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N  is  termed  the  volume  fraction  index  which  provides  the  material  variation  profile  through  the 
plate  wall  thickness,  (0  <  <  qo)  .  A  pictorial  representation  of  the  distribution  of  the  constituent 

materials  are  shown  below  in  Fig.  2. 

Case  (2):  The  phases  vary  non-symmetrically  through  the  wall  thickness,  and  in  this  case  there  is 
full  ceramic  at  the  outer  surface  of  the  plate  wall  and  full  metal  at  its  inner  surface.  For  this  case, 
(z)  can  be  expressed  as 


Vc(^)  = 


^  h  + 


(6) 


Below  is  a  pictorial  representation  of  the  antisymmetric  case  shown  in  Fig.  3. 

It  should  be  noted  that  in  contrast  to  case  (2),  where  there  exists  coupling  between  stretching 
and  bending,  such  coupling  is  not  present  for  the  symmetric  case  (1).  Also,  for  the  purposes  of 
simplicity  the  Poisson’s  ratio  will  be  assumed  to  be  constant  throughout  the  plate  structure.  From 
Eqs.  (l)-(6),  the  effective  material  properties  of  a  FG  plate  can  be  expressed  for  the 
asymmetric  case  as 


[E{z),  p{z)]  =  {E 

cm  ’  Pcm  ) 


h  +  2z 

2h 


Pm) 


(V) 


and  for  the  symmetric  case  as 


[E(z),/j(z)  ’  Pcm  ) 


A  _  \Nc 


vV2y 


l  +  sgn(z) 


,A/2, 


N , 


1  -  sgn(z) 


v(z)  =  V  , 


^^m  ’  Pm  )  (^) 


(9) 


Where 


P'cm  ~  P'c  P'mi  Pcm  ~  Pc  Pm- 


(10) 
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3.  KINEMATIC  EQUATIONS 

3.1  The  3-D  Displacement  Field 

Consistent  with  the  classical  plate  theory  [6],  the  distribution  of  the  3-D  displacement 
quantities  through  the  wall  thickness  can  be  expressed  as 

u=uq-zwq^^,  v  =  VQ-zwQ^y,  w^Wq.  (lla-c) 

Within  these  equations,  (u,v,w)  are  the  3-D  displacement  quantities  along  the  {x,y,z) 
directions,  respectively.  While,  uq,Vq,  Wq  are  the  2-D  displacement  quantities  of  the  points  on 
the  mid-surface. 

3.2  Non-Linear  Strain-Displacement  Relationships 

The  nonlinear  strain  displacement  relationships  across  the  plate  thickness  at  a  distance  from  the 
mid-surface  take  the  form  [2,  3,  6] 


'XX 


^yy  ^,y  ^ ^  '^^xy  fxp 


=  M  „  d-U  r  +  W  rW  v  ■ 
r>y  9-^  9-^  '>y ' 


(12) 


-■XZ 


0. 


Substitution  of  Equations  (lla-c)  into  Equations  (12)  results  in  the  strain  measures  across  the 
plate  thickness  in  terms  of  the  2-D  displacement  quantities  of  the  mid-surface  of  the  plate 
expressed  as 


r  1 

L(o)l 

L(i)l 

^xx 

^yy 

>  —  < 

>  +  z< 

4' 

Jxy  ^ 

I4J 

(13) 


Where 


{e 


(0) 

XX 


(14a) 
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i  ^0,xx  ^0,}y  ^^0,xy}-  (14b) 

In  the  above  expressions,  are  referred  to  as  the  membrane  strains  and 

referred  to  as  the  flexural  bending  strains  which  are  also  known  as  the 

curvatures. 

4.  CONSTITUTIVE  EQUATIONS 

The  stress-strain  relationships  for  a  state  of  plane  stress  is  expressed  as  [10] 


^  XX 

'Qn 

012 

0  ■ 

^xx 

^yy 

>  = 

Qn 

022 

0 

^yy 

^xy 

0 

0 

066  _ 

Txy  ^ 

^xz  ^yz  ^  zz  (!•  (15) 


The  material  stiffnesses,  Qy{z),{i  =  1,2,6)  are  given  by  [9,  10] 

011=022  =^,212=^,266=:^.  016=026=0. 
1-v^  l-v^  2[1  +  V] 

The  standard  force  and  moment  resultants  of  a  plate  are  defined  as 

•A/2 


(16) 


rniz 

(Nij  ’ )  "  J_  ’  iiJ  = 


(17) 


With  the  use  of  Equations  (12)-(17),  the  stress  resultants  and  stress  couples  are  related  to  the 
strains  by  [3] 


fwi 

'[A]  [B] 

1{M}J 

[B]  [Dl 

W"}J 

(18) 


in  which  [A],  [B],  and  [D]  are  the  respective  in-surface,  bending-stretching  coupling,  and 
bending  stiffnesses.  For  the  case  of  symmetric  EG  Plates,  [B]=0,  since  there  is  no  bending- 
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Stretching  coupling.  The  global  stiffness  quantities,  5^,  and  ,  (/,  7=1,2,6)  are  defined 
respectively  as 

(Ay.Bij.Dij)  =  a,-(l,z,zb<i'z .  (19) 

In  view  of  Equations  (7),  (8),  (9),  (16),  and  (19),  the  global  stiffness  quantities  are  determined 
as 

[(4d^22)»(^1D^22X(Ad^22)]  = - A(^b^2»^3), 

1  -  V 

(^12  ’  ^12  ’  A2  )  = - A  (^1  ’  ^2  ’  ^3  ) .  (20a-c) 

1-v^ 

(^66’^66’^66)  - r(^D^2’^3)’ 

2(1  +  v) 

where  for  the  asymmetric  case, 

Ei=E^h[l  +  l/(N  +  l)],  E2^E^h^[l/(N  +  l)-l/(2N  +  l)]  (21) 

£■3  =  / 12  +  /(TV  +  3)  - 1 /(TV  +  2)  + 1 /(4iV  +  4)] 

And  for  the  symmetric  case, 

El  =E^h[\  +  \l{N  +  l)'\,  E2=0,  E2  ^  E^h^[l/l2  +  l/  4(N  +  3)].  (22) 

5.  GOVERNING  EQUATIONS 

Hamilton’s  principle  is  used  to  derive  the  equations  of  motion  and  the  boundary  conditions.  It 
is  formulated  as 

SJ^d{‘\K-U  +  V)dt^O,  (23) 

Jto 
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where  to  and  t^are  two  arbitrary  instants  of  time.  U  denotes  the  strain  energy,  V  denotes  the  work 
done  by  surfaee  traetions,  edge  loads,  body  forees,  and  viseous  damping  forees,  and  K  denotes 
the  kinetie  energy  of  the  3-D  body  of  the  strueture,  while  5  is  the  variational  operator.  In 
Equation  (23),  the  strain  energy  is  given  by 


dU^ 


2 


xx^^xx  ^  yy^^yy  ^  xy^y xy^^^^^Q 


yy'^'^yy 


xy'^rxy . 


(24) 


where  Qq  denotes  the  mid-surfaee  area  of  the  panel.  The  work  done  by  external  loads  is 
expressed  as 


dV 


n 

=  I*  [Pf{x,y,t)-2nw^wdQ.Q+[ 


,  a,„,Svdzdx  + 

h  yy 


h 
2 
y- 


II 


^  a  j^dudzdy. 


(25) 


In  the  above  expression  Pf(x,y,t)  is  the  transverse  pressure,  jj.  is  the  viseous  damping 


eoefficient  per  unit  area  of  the  plate,  (cr^^^ ,  cr^ )  are  the  specified  stress  components  along  the 

plate  edges,  and  (du,  dv)  are  the  virtual  displacements  in  the  normal  and  tangential  directions, 
respectively,  along  the  plate  edges.  Considering  only  the  transversal  inertia  of  the  structure,  the 
variation  in  kinetic  energy  is  expressed  as 

dK  -  I*  iQwdwdPlQ,  (26) 

JQq 


where 


2 


(27) 


and  p{z)  is  the  mass  per  unit  volume. 
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Considering  Equation  (23)  in  conjunction  with  Equations  (24)-(27)  ,  along  with  the  constitutive 
equations  (15),  the  strain-displacement  relationships,  equations  (14),  carrying  out  the  integration 
through  the  thickness,  integrating  by  parts  whenever  feasible,  using  the  expression  of  global 
stress  resultants  and  stress  couples,  while  retaining  only  the  transversal  load,  transverse  inertia, 
and  transverse  damping,  and  invoking  the  arbitrary  and  independent  character  of  variations 
duQ,dvQ,dwQ,dwQ^^,dwQy,  one  obtains  the  equations  of  motion  and  as  a  by-product  the 


boundary  terms  or  conditions.  This  results  in  three  equations  of  motion  in  terms  of  the  global 
stress  resultants  and  stress  couples  and  four  boundary  conditions  along  the  plate  edges.  These 
equations  of  motion  and  boundary  conditions  can  be  respectively  expressed  as 

^0-  ^xx,x+^xy,y  (28) 


8v)q  :  V  -F  N- 


Swn  : 


yy,y  "'xr,x 


0 


^xx,xx  +^^xy,xy  '^^yy,yy  ^ xy^0,y^ X  yy^0,y 

^xy^Q,x)y  +  “  2^  =  IqWq 


(29) 


(30) 


The  boundary  conditions  become 
Along  the  edges  x  =  Q,  L\ 

=0 

^xx,x  ~^'^^xy,yy  ~^^xx^0,x  ^ xy^0,y 

Along  the  edges  y  =  0,  L2 

Ny,=^ 


or 

0 

II 

0 

s 

(31) 

or 

c 

0 

II 

0 

(32) 

or 

Wo,x  =  0 

(33) 

0 

or 

0 

II 

0 

(34) 

or  wq  =  0  (35) 
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o 

II 

or 

^>0 

=  0 

(36) 

Myy  0 

or 

^0,y 

=  0 

(37) 

^)y,y  ^^xy,xx  +  ^xy^Q,x  +  ^yy^Q,y  = 

or 

Wo 

=  0 

(38) 

For  the  case  of  all  edges  simply  supported  and  freely  movable,  the  boundary  conditions  are  as 
follows: 

wo=M^=N^=Q,N^=Nl,onx  =  0,L^  ^9^ 

wo  =Myy  =Ny^  =  ^ ,  N  yy  =N*yy  On  ^  =  0,  ^2  (40) 

It  should  be  mentioned  for  clarification  sake  that  for  compressive  edge  loading, 

nIc  =-Nl,,m,dN*yy  =-N%,. 

AIR-BLAST  LOADING 

With  the  ever  increasing  demands  for  increased  safety  for  the  warfighter  in  the  field  to  operate 
structurally  sound  vehicles  in  the  event  of  an  lED  or  some  other  type  of  explosive,  it  is 
imperative  that  an  understanding  of  the  structural  response  of  various  components  within 
military  combat  vehicles  under  an  explosive  blast  be  understood  so  that  measures  can  be  taken 
from  a  design  standpoint  to  ensure  the  durability  and  survivability  of  these  components.  To  begin 
to  achieve  this  understanding,  the  type  of  explosive  loading  considered  here  is  a  free  in-air 
spherical  air  burst.  Such  an  explosion  creates  a  spherical  shock  wave  which  travels  radially 
outward  in  all  directions  with  diminishing  velocity.  The  form  of  the  incident  blast  wave  from  a 
spherical  charge  is  shown  in  Fig.  4.  Where  Psq  is  the  peak  overpressure  above  ambient  pressure, 
Pq  is  the  ambient  pressure,  is  the  time  of  arrival,  tp  is  the  positive  phase  duration  of  the  blast 
wave,  and  t  is  the  time.  The  waveform  shown  in  Fig.  4.  is  given  by  an  expression  known  as  the 
Friedlander  equation  and  is  given  as 
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Pt  (0  =  {PsQ -PQ)W-{t-ta)ltp]  [-a(f  -t  a)  It  pi  (41) 

where 

P50  =1772/Z^-114/Z^+108/Z.  (42) 

1/3 

In  equations  (41)  and  (42),  Z  is  known  as  the  sealed  distanee  given  by  Z  =  RIW  with  R 
being  the  standoff  distanee  in  meters  and  W  being  the  equivalent  eharge  weight  of  TNT  in  terms 
of  kilograms.  Also,  a  is  known  as  the  deeay  parameter  whieh  is  determined  by  adjustment  to  a 
pressure  curve  from  a  blast  test. 

For  the  conditions  of  standard  temperature  and  pressure  (STP)  at  sea  level,  the  time  of  arrival 
and  the  positive  phase  duration  tp  can  be  determined  from  [5] 

tltx=RIRx={WIWxf'^.  (43) 

ti  represents  either  the  arrival  time  or  positive  phase  duration  for  a  reference  explosion  of  charge 

weight  W\ ,  and  t  represents  either  the  arrival  time  or  positive  phase  duration  for  any  explosion 
of  charge  weight  W .  The  determination  of  the  standoff  distance  for  any  charge  weight  W  follows 
a  similar  reasoning.  The  application  of  these  relationships  is  known  as  cube  root  scaling.  It 
should  be  understood  that  in  applying  these  relationships  that  the  standoff  distances  are 
themselves  scaled  according  to  the  cube  root  law. 

6.  SOLUTION  METHODOLOGY 

To  satisfy  the  first  two  equations  of  motion,  equations  (28,29),  a  stress  potential  will  be  utilized 
which  allows  the  in-plane  stress  resultants  to  be  expressed  by  letting 

■l^xx  ~  V,yy  5  yy  ~  V,xx  ’  ^^xy  ~  ~V,xy  •  (44) 
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The  third  equation  of  motion,  equation  (30),  can  be  expressed  in  terms  of  two  unknown 
variables,  the  stress  potential  ^and  the  transverse  displacement  Wg.  To  accomplish  this,  a  partial 
inversion,  of  equation  (18),  the  details  of  which  are  not  presented  here,  needs  to  be  carried  out. 
Performing  a  partial  inversion  results  in  [4] 

U'U}  J 

where 

[A*]  =  [A]-\  [B*]  =  ^AT\b]  ,  -[S*f  =[S][^rl 

[D  ]  =  [D]-[B][A\-\bI  (46) 

Using  Equations  (44),  (45),  and  (46)  and  simplifying,  equation  (30)  takes  the  form 

-{(p^xx^Q^yy  -^(P,xy^Q,xy  +  (P,yy^Q,xx)  + +1^Wq  =  Pf  (47a) 

where 


[A*] 

-[B*f 


[B*] 

[D*] 


{N} 


1451 


D  = 


E1E3  -Ej 

E,{\-v^) 


e  =  -^. 


(47b,c) 


This  gives  one  governing  equation  with  two  unknowns,  wq  and  (p .  One  more  equation  is  needed 
in  terms  of  the  same  unknowns  which  will  give  two  equations  in  terms  of  two  unknowns  which 
can  then  be  solved.  This  will  come  from  the  compatibility  equation.  By  eliminating  the  in-plane 
displacements  from  the  strain-displacement  relationships,  equations  (14)  the  relationship 
between  the  in-plane  strains  and  the  transversal  deflection  known  as  the  compatibility  equation 
can  be  shown  to  be  given  by 


0  ,0  0  _  2 

^xx,yy  yxy,xy  ^0,xy  ^0,xx'^0,yy 


(48) 


In  view  of  equations  (44)  ,  (45),  and  (46),  the  compatibility  equation  is  expressed  as 
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^^(p  =  Ex{wQ^:^-WQ^^WQ^yy).  (49) 

In  equations  (47a)  and  (49),  +5^5^^  whereV  is  referred  to  as  the  Laplaeian 

operator. 

Equations  (47a)  and  (49)  are  the  basie  governing  equations  used  to  investigate  the  struetural 
response  of  EG  plates  under  external  exeitation  loading.  Eor  the  purposes  of  this  paper,  from  this 
point  forward,  the  thermal  terms  will  be  disearded.  To  this  end,  to  solve  equations  (47a)  and  (49), 
the  approaeh  adopted  from  [2]  will  be  utilized.  In  this  respeet,  the  following  funetional  forms  are 
assumed  for  Wq  and  ^[2]. 


wq  (x,  j,  f (0  sin  sin  ju„y 
~  ^mn  (^)cos  2/l^x  +  (t^ COS  +  Cjfifi  (t)eos  2/l^xeos  'T.fiy^y  + 

/  N  .  .  1  *  2  1  *  2  (50a,b) 

vvsin ^sin jF  +  —  jF  +  — , 


where  -rmi:ILi,id^  -nnl L2,m,n-\,2,...  are  the  number  of  half  waves  in  the  xand  y 
direetions,  respeetively,  and  is  amplitude  of  defleetion.  Also, 

^mn  (0^  (0^  Cwu  (0^  (0  ^^e  ooeffioients  to  be  determined.  By  substituting  equations 

(50a,b)  into  equation  (49),  the  ooeffioients  A^„(i),  5^„(i),  C^„(i),  and  Z)^„(i)  are  determined 
as 

^mn  (0  =  El  w^„  (t)nl  /  322^ ,  it)  =  El  (02^  /  32//^ ,  (t)  =  (0  =  0-  (51) 

As  a  result  of  the  assumed  form  for  the  Airy’s  potential  funotion,  (p  the  following  relationships 
expressed  mathematioally  as 


C^l 

1  * 

ch 

1  * 

Jo 

II 

O 

> 

II 

K) 

1  V',xx 

Jo 

II 

o 

to 

II 

(52) 
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implies  that  Ni  j ,  N22  acquire  the  meaning  of  average  in-plane  compressive  edge  loads. 
Moreover,  by  assuming 

Pt{x,y  ,t)  —  f^ny  ■>  (53a) 

which  implies  through  integration  of  both  sides  over  the  plate  area  that 

4  cPi  f  A 

= -n- \  Pt{x,y,t)?,mX^x?,mii„ydxdy .  (53b) 

L1L2  Jo  Jo 


Letting, 

Pt (x, y, t)  =  Pt (0  =  {Pso  - Po )[1  -{t-ta)ltp]Qy^ [-a{t -ta)ltp]  (53c) 


and  integrating  gives 

Pmn  (0  =  (0  /  •  (m,  n)  =  (1, 1)  (53d) 

Introduction  of  equations  (50a,b)  and  (51)  into  equation  (47a)  and  retaining  the  resulting 
equation  along  with  the  unsatisfied  boundary  conditions  in  the  energy  functional  and  applying 
the  Extended  Galerkin  technique  results  in  the  following  ordinary  nonlinear  differential  equation 
governing  the  structural  response  of  EG  plates,  under  external  excitation  which  is  expressed  as 

•  2  3  ^ — 

^mn  (0  +  2A 

mn  ^mn  ^mn  (0  "b  ^  mn  ^mn  (0  "b  ^mn  (O^mn  (0  Pmn  it),  (54) 

where  WynniP)  represents  the  amplitude  of  deflection  of  the  plate  as  a  function  of  time, 

Pmn  it)  =  167)  it)  /  IqTI^  ,  (O^n  =  /  Iq  is  the  natural  frequency  of  the  EG  plate. 


^mn  —  di ^t)^mn  1^  the  non-dimensional  damping  factor,  and  -  ^l(^w +/^«)/167o  •  11 
should  be  noted  that  at  the  center  of  the  plate  {x,y)  =  /  2,  L2  !  2) ,  Wjyj„it)  is  equal  to  the 
maximum  deflection  of  the  plate.  In  these  latter  expressions, 

Kmn  =mi  +My  +N*^2i  +N*yy^l  .  (55) 
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Equation  (54)  is  a  nonlinear  equation  in  terms  of  the  plate  deflections  as  a  function  of  time.  It  is 
interesting  to  note  that  equation  (54)  is  very  similiar  to  Buffing’s  Equation.  To  obtain  the  plate 
deflections  as  a  function  of  time,  equation  (54)  is  solved  using  the  Eourth-Order  Runge-Kutta 
Method  with  zero  initial  conditions. 

6.1  Special  Case  -  Constant  Uniform  distributed  Transverse  Pressure 

Eor  validation  purposes  later  on,  the  special  case  of  a  constant  uniformly  distributed  transverse 
pressure  acting  on  the  top  surface  of  the  plate,  independent  of  time,  for  a  square  plate  will  be 
shown  [1].  Eor  a  square  plate  under  a  uniform  pressure  independent  of  time,  w^„(t)  =  0,  and 

Wmn  (0  =  0  •  As  a  result,  equation  (54)  reduces  to 


2  3 

^mn  ^mn  ^mn  ^mn  ~  ^mn  • 


(56) 


Neglecting  tensile  or  compressive  loading,  equation  (56)  can  be  written  as 


'^\  4  4  3  2  2  2 

,  ^  f^n  )  ^mn  ~  ^mn 

16 


Where, 


(57) 


P  = 

mn 


\6P 


(58) 


mnn 


Equation  (57)  is  a  nonlinear  algebraic  equation  which  can  be  solved  by  a  typical  numerical 
technique  for  the  load-deflection  relationship.  This  will  be  shown  graphically  in  the  results 
section. 

7.  ELASTIC  STRESS  ANALYSIS 

With  the  transversal  deflections  in  hand,  as  a  function  of  time,  and  with  the  use  of  equations  (13), 
(14a,b),  (15),  and  (16),  the  constitutive  equations  become 


E(z) 


1-v" 


W0,x  +  tr  W0,x  +  V  +  tr  W0,y  -  +  ^0,yy  ) 

z  y  z  j 
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Gyy  (X,y,Z,t)  = 


E(z) 


a^{x,y,z,t)  = 


1-v^ 

E{z) 


1-  -2 
2 


V(),y+7:^Q,y+V 


^  2 
^0,x  +  T 

V  2 


2(1 +  v) 


-Z{WQ^yy  +VW0,xr) 

U0,y  +  +  ^0,x^0,y  “  _ . 


(59a-c) 


At  this  point,  all  thermal  terms  have  been  discarded  since  the  thermal  response  is  not  the 
subject  of  this  paper.  In  order  to  determine  the  stresses,  the  expressions  for  Uq  and  uq  in  terms 
of  the  transversal  amplitude  of  deflection,  (t)  needs  to  be  determined.  These  relationships 
can  be  determined  from  the  global  constitutive  equations  (18).  Expanding  equation  (18)  gives, 

^xx  +^12^jy^  +^11^«  +^12^>y 


^yy  -  +^22^j^^  +^12^ix  +^22^jy  • 


(60a,b) 


Solving  equations  (60a,b)  simultaneously  for  ,  and  in  terms  of  ,  Ny^ ,  ,  and 


and  utilizing  equations  (14)  and  (44)  gives. 


1 


1 


^0,x  ^  ^^,yy  ^V,xx^~^  ^  ^0,xr 


1  .  ^  ^2  1  2 
^0,y  -  iV,xx  ^V,yy^^  ^0,yy  2^®’-^’ 


(61a,b) 


With  the  use  of  equations  (20a-c),  (50a,b),  and  (51),  and  integrating  results  in. 


Uq 


-  Sin  xi 


8 


-X-- 


16A„ 


22™  .  2 


sm  Uny-v 


V 


^  (0^2 2^  2^xsin/i„y 

E\ 


+  ~(^XC  -VEfy^)X  +  Cl 

E\ 


Page  16  of  37 

UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


_  '^mn  (^)Mn  ^mn  •  q  -  2  i  ^mn  i  v 

Vq  ^  1  z'  9  V  +  SUl /l^XCOS jp 

8  l6//„  & 

r-n  V  y  i 


*  * 


+  — -vN^)y  +  C2 

^1 


(62a, b) 

Finally,  substitution  of  equation  (50a),  and  equations  (62a,b)  into  equations  (59a-o)  gives  the 
expressions  of  the  stresses  as. 


,  ,  E{z)  2.9  0  (  ^2^.2  2\  .9  •  , 

= - yi - z - w^n(0cos2^„y-  z-—  +vfi„  ^^„(t)smX^xsm^„y 

l-v^  [  8  ^  ' 

+nJ 


Oyy{X,y,Z,t)=- 


Eiz^  (^1  ^  ^  Et  ^ 

y  I  Z  (0  —  Z  \^n  j'^wn  (0 


El  '  Ej  ' 

(■^’ F’ -^5  0  ~  ^  (0 

1  +  v  K  fij  j 


(63a-e) 


Where, 


T7  ’  yy  r 


(64a,b) 


are  the  nondimensional  eompressive  edge  loadings. 

7.1  Special  Case  -  Constant  Uniformly  Distributed  Transverse  Pressure 

Under  the  same  assumptions  used  in  seetion  6.1,  the  in  plane  axial  stress,  ,  is  given  as 


^xx 


E{z)  \( 2  f  ^2\,2  29 

o  I  o  ^mn  +  „  (■^m  )^mn  f  • 

1-v^  8  I  ^iJ 
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It  should  be  mentioned  that  for  a  square  plate,  =  ayy  .  Also,  =  0  at  the  center  of  the 

plate.  Once  the  deflections  are  determined  from  equation  (57),  they  are  then  substituted  into 
equation  (65)  to  determine  the  load-stress  relationship.  This  will  be  shown  graphically  in  the 
results  section. 

8.  RESULTS  AND  DISCUSSION 

To  validate  the  present  theory,  comparisons  are  made  with  Woo  and  Meguid  [1]  who 
considered  a  transverse  uniform  pressure,  independent  of  time,  acting  on  the  surface  of 
asymmetric  functionally  graded  (FG)  plates  and  shallow  shells.  Figures  5  and  6  show 
comparisons  between  the  present  results  and  Woo  and  Meguid.  The  geometrical  and  material 
properties  used  for  the  comparison  are  given  as 

Zj=0.2w,  A  =  0.01m,  y/  =  l 
Aluminum :  =  70  GPa,  v  =  0.3 

Alumina  :  =  380  GPa,  v  -  0.3 . 

Figure  5  reflects  the  relationship  between  the  dimensionless  central  deflection  and  the  load 
parameter  Q  =  PfLj  j .  Inspection  of  both  results  (Present  and  Woo  and  Meguid)  reveals 
somewhat  of  a  marginal  difference.  This  marginal  difference  seems  to  increase  as  the  load 
parameter  Q  becomes  larger.  Also  the  differences  appear  larger  for  pure  metal  (Aluminum)  as 
opposed  to  pure  ceramic  (Alumina).  In  both  cases,  with  larger  deflections,  comes  larger 
differences.  The  differences,  most  likely  come  from  the  fact  that  Woo  and  Meguid  assumed  a 
different  solution  process  which  consisted  of  a  combination  of  series  solutions  and  or  Fourier 
series  representations,  for  their  deflection  profile,  in  combination  with  numerical  methods.  For 
the  present  solution  process,  a  one  term  Navier  type  solution  was  utilized  for  the  deflection 
profile.  In  addition,  the  Galerkin  Method  which  is  an  approximate  method  was  used  to  facilitate 

Page  18  of  37 

UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


the  governing  equation  with  was  then  solved  using  the  Runge-Kutta  mathod.  Two  different 
solution  processes  were  utilized  between  the  present  author  and  Woo  and  Meguid. 

2  /  2 

Figure  6  highlights  the  relationship  between  the  dimensionless  axial  stress  o  =  OxL\  j  E^h 
and  the  dimensionless  thickness  4’  =  z//z  for  a  fixed  value  for  the  load  parameter  Q  =  -400. 
Since  there  are  differences  existing  with  the  deflection-load  history,  their  will  be  differences 
revealed  with  the  stress  distribution.  By  observation,  it  can  be  seen  that  some  marginal 
differences  exist.  Again,  this  is  attributed  to  the  fact  that  two  unrelated  approaches  were  utilized 
to  arrive  at  the  same  final  results. 

To  illustrate  the  present  approach,  a  ceramic-metal  functionally  graded  plate  consisting  of  Ti- 
6A1-4V  and  Aluminum  Oxide,  with  the  following  material  properties,  were  considered  for  the 
numerical  results  presented. 

Ec  =  320.24  GPa,  =  3750  kg  /  =  0.26 

^^=105.7GPa,  =4429kg/m^  =0.2981 

^ave  =  0-2791 

The  geometrical  properties  used  for  the  FG  Plate  are  Tj  =lm,  if/  =  Li/L2  =1,  and  unless 
otherwise  stated  h  =  0.0254  m,  the  decay  factor,  or  =  0.5,  and  the  halfwaves,  {m,n)  =  (1,1).  In 
addition,  the  following  reference  values,  in  Table  1,  were  utilized  to  determine  the  time  of  arrival 
and  positive  phase  duration  [5]. 

Figure  7  depicts  the  central  2-D  stress-time  history  of  an  asymmetric  FG  plate,  on  the  top 
surface  of  the  plate,  with  damping  and  an  aspect  ratio  of  1.25.  both  stress-time  histories  have  the 
same  frequency  but  different  amplitudes.  It  is  revealed  that  for  aspect  ratios  greater  than  one  a yy 

has  larger  amplitudes.  Although  not  shown  here,  when  the  aspect  ratio  is  equal  to  1,  the  stress- 
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time  history  for  ,  Oyy  are  the  same  (same  frequeney  and  amplitudes).  When  the  aspeet  ratio  is 
less  than  1 ,  the  amplitudes  of  the  stress-time  history  of  is  larger  than  for  Oyy  with  identieal 
frequencies.  In  Figure  8,  the  central  stress-time  histories  of  the  stress  component,  a ^  on  the  top 

surface  of  the  plate,  for  various  volume  fraction  indexes,  is  depicted.  Three  comparisons  are 
made  between  fully  ceramic,  fully  metal,  and  FG  with  ceramic  and  metal.  As  the  plate 
construction  transforms  from  fully  ceramic  to  fully  metal,  it  does  so  with  decreasing  frequency 
and  increased  stress  amplitude.  In  other  words,  the  fully  ceramic  case  has  the  highest  frequency 
due  to  its  increased  stiffness  but  with  a  lower  stress  amplitude.  In  contrast  the  fully  metal  case 
has  the  highest  stress  amplitude  and  the  lowest  frequency  due  to  its  decreased  stiffness  when 
compared  with  ceramic.  This  figure  shows  that  for  some  applications  ceramic  may  be  more 
beneficial  to  adopt.  Figure  9  depicts  the  damped  stress-time  histories,  of  an  asymmetric  plate, 
for  the  2-D  stress  component,  <7^  at  three  different  central  locations  which  include  the  top 

surface,  the  middle  surface,  and  the  bottom  surface.  It  can  be  seen  that  on  the  central  top  surface 
(ceramic)  the  stress-time  history  starts  out  in  compression  along  with  the  top  surface  in  tension. 
It  is  also  revealed  that  the  central  stress  amplitude  on  the  top  surface  (ceramic)  is  larger  than  for 
the  bottom  surface  (metal).  The  ceramic  bears  a  large  portion  of  the  stress  since  the  ceramic  acts 
as  a  restraint  in  preventing  the  metal  surface  from  a  large  amount  of  deformation  and  or  stress 
build-up.  The  central  stress-time  history  for  the  middle  surface  is  very  small  in  amplitude, 
although  not  zero  (not  a  neutral  surface).  Figure  10  is  the  counterpart  of  figure  9  for  the  case  of  a 
symmetric  FG  plate  with  ceramic  on  the  top  and  the  bottom  surface  with  metal  in  between.  It  can 
be  seen  that  the  central  stress-time  histories  for  both  the  top  and  bottom  surfaces  are  almost 
reflections  of  one  another.  Both  have  the  same  frequency  but  are  out  of  phase  with  each  other. 
Although  almost  symmetrical,  it  should  be  noted  that  the  stress  on  the  middle  surface  is  not  zero 
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in  contrast  to  beam  theory,  the  middle  surface  is  not  a  neutral  surface  but  very  close.  The  stress¬ 
time  history  is  negligible  when  compared  to  the  top  and  bottom  surfaces.  Figure  1 1  compares  the 
central  stress-time  history  for  on  the  bottom  surface  for  both  the  asymmetric  and  symmetric 

FG  constructions.  While  both  histories  exhibit  the  similar  behavior,  the  symmetrie  ease  generates 
larger  amplitudes  when  it  is  in  both  eompression  and  tension.  The  plausible  explanation  for  this 
behavior  is  that  both  ceramie  surfaees  whieh  by  nature  are  stiffer  than  the  middle  surface  being 
metal  act  as  a  constraint  preventing  the  metal  from  an  attributable  amount  of  deformation.  The 
eeramie  bears  most  of  the  stress  with  little  deformation.  Figure  12  is  the  eounterpart  of  figure  1 1 
comparing  the  symmetries  at  the  middle  surface.  Again  although  the  stress  profiles  are  small 
they  are  not  zero.  The  middle  surfaee  is  not  a  neutral  surface.  Finally,  in  figure  13,  whieh  is  the 
counterpart  of  both  figures  1 1  and  12  at  the  top  surface  of  the  plate,  it  is  shown  that  there  is  only 
a  negligible  difference  in  the  central  stress-time  histories  between  the  asymmteric  and  symmetric 
case.  It  appears  that  in  both  cases  the  ceramic  bears  the  same  amount  of  stress  irrespective  of  the 
symmetry. 

As  a  final  point  of  discussion,  a  comment  or  two  should  be  made  about  permanent 
deformation  and  fracture.  The  theory  presented  here  is  based  on  elasticity  theory  whieh  does  not 
take  plastie  deformation  or  fraeture  into  aecount.  It  is  assumed  that  the  behavior  of  the  plate 
remains  within  the  elastic  limit.  To  determine  if  the  elastic  limit  is  exceeded  by  way  of 
deformation  or  induced  stresses,  appropriate  verification  and  validation  (V&V)  would  need  to  be 
carried  out  through  laboratory  testing  per  the  required  material  phase  (ceramic  and  metal) 
distribution  throughout  the  structure.  Another  option  ineludes  (V&V)  against  finite  element 
modeling  and  simulation  (M&S)  software  whieh  may  inelude  a  more  accurate  material  model 


Page  21  of  37 

UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


such  as  the  Johnson-Cook  material  model  based  on  a  high  strain-rate  of  behavior  of  loading 
whieh  is  not  captured  in  the  elasticity  theory. 

9.  CONCLUDING  REMARKS 

A  rigorous  treatment  of  functionally  graded  plates  with  grading  in  the  transverse  direetion  for 
two  different  types  of  funetionally  graded  symmetries  (asymmetrie  and  symmetrie)  has  been 
studied.  Validations  with  a  simpler  case  eonsisting  of  a  transverse  uniform  load  found  in  Woo 
and  Meguid  [1],  has  been  made.  It  has  been  shown  that  damping  has  an  important  effeet  when  it 
eomes  to  the  attenuation  of  the  deflections.  It  has  also  been  shown  that  other  faetors  sueh  as  the 
magnitude  of  the  volume  fraetion  index,  the  aspeet  ratio,  the  loeation  of  the  stress,  and  most 
importantly,  the  symmetry  of  the  transverse  grading  throughout  the  structure  plays  an  important 
role  in  the  eentral  stress-time  history  of  the  strueture. 

It  has  also  been  shown  that  the  eeramie  carries  a  large  part  of  the  axial  stress  while  aeting  as  a 
eonstraint  in  the  deformation  of  the  metallie  phase.  Also,  the  neutral  surfaee  does  not  reside  in 
the  middle  surface  of  the  plate.  In  addition,  it  is  apparent  that  the  stress-time  history  of  the 
metallie  phase  exhibits  larger  amplitudes  of  stress  but  lower  frequeneies  in  eontrast  to  eeramie 
whieh  exhibits  the  opposite  behavior  due  to  its  stiffness. 

From  a  design  standpoint,  depending  on  the  design  requirements,  it  would  be  appropriate  to 
state  that  integration  of  FG  materials  in  plate-type  structures  would  benefit  the  structure  response 
of  the  strueture.  It  is  hoped  and  realized  that  this  present  study  presented  here  will  give  insight 
into  some  of  the  factors  that  can  play  an  important  role  in  the  struetural  response  of  funetionally 
graded  plates  and  fill  in  some  of  the  fundamental  missing  gaps  within  this  subjeet  area. 
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FIGURE  CAPTIONS 


Fig  1 .  A  simply  supported  functionally  graded  plate  shown  in  2-D  under  an  explosive  blast. 

Fig  2.  Distribution  of  the  constituent  materials  through  the  plate  thickness  for  the  symmetric 
case. 

Fig  3.  Distribution  of  the  constituent  materials  through  the  plate  thickness  for  the  antisymmetric 
case. 

Fig  4.  Incident  pressure  profile  of  a  blast  wave. 

Fig  5.  Validation  of  the  dimensionless  central  deflection  versus  uniform  transverse  load. 

Fig  6.  Validation  of  the  dimensionless  axial  stress  along  the  thickness  z  at  the  center  of  the  plate 
under  the  load  parameter  Q  =  -400 

Fig  7.  2D  stress  components  versus  time  for  an  asymmetric  FG  plate  (V=0.5). 

Fig  8.  The  axial  stress  component  along  the  x  direction  versus  time  for  various  values  of  the 
volume  fraction  index. 

Fig  9.  The  axial  stress  component  along  the  x  direction  versus  time  at  various  locations 

throughout  the  thickness  at  the  center  of  the  plate  for  an  asymmetric  FG  Plate  (V=0.5). 

Fig  10.  The  counterpart  of  figure  9  for  a  symmetric  FG  plate  (V=0.5). 

Fig  1 1 .  The  axial  stress  component  along  the  x  direction  versus  time  for  both  the  asymmetric 
and  symmetric  FG  plate  on  the  bottom  surface  at  the  center  of  the  plate  (V=0.5). 

Fig  12.  The  axial  stress  component  along  the  x  direction  versus  time  for  both  the  asymmetric 
and  symmetric  FG  plate  on  the  middle  surface  at  the  center  of  the  plate  (V=0.5). 

Fig  13.  The  axial  stress  component  along  the  x  direction  versus  time  for  both  the  asymmetric 
and  symmetric  FG  plate  on  the  top  surface  at  the  center  of  the  plate  (V=0.5). 
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TABLES 


Table  1,  Airblast  Parameters  Versus  Distance  for  a  One  Kilogram  (iFi )  TNT  Spherical 
Air  Burst  [5], 


Standoff  Distance,  i?i(m) 

Arrival  Time,  (msec) 

Positive  Phase  Duration,  t (msec) 

1.0 

0.532 

1.79 
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FIGURES  (In  Numerical  order) 
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Dimensionless  Thickness, 
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